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Abstract 

We investigate the impact of finite temperature effects on axions in the context of cos- 
mology. The temperature dependence of the decay constant is modeled analogous to pions. 
For the two interesting cases considered here, we find that the temperature effects do lead 
to changes relevant for detailed and precise abundance and rate calculations. We also find 
that the axion decoupling temperature starts showing large deviations for larger values of the 
axion decay constant. 
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1 Introduction 



Axions, the pseudo Goldstone bosons, associated with the Peccei-Quinn symmetry, U(1)pq, remain 
the most popular solution to the strong CP problem P3 |21 El H|- Axions, QCD type or some 
other variant, appear in many extensions of the standard model. The mass and coupling of the 
axion is inversely proportional to the scale at which U(1)pq is broken, denoted by /„, and here 
after referred to as Peccei-Quinn (PQ) scale or axion decay constant interchangeably. This has 
important implications for axion cosmology because axions can only be produced thermally when 
the temperature of the universe drops below f a . The PQ scale, f a , is constrained by astrophysics, 
cosmology and laboratory experiments. For a review of latest bounds on f a refer to E] , while 
a more comprehensive review/collection of (earlier) bounds on f a can be found in |7| IH} 19*1 ITUl ITT]. 
Further, axions are also one of the theoretically well motivated contenders for the cold dark matter 
in the universe. The DAMA collaboration has also explored the possibility of scalar and pseudo- 
scalar interpretation of their results ^2]- More recently, results from the PVLAS experiment 
indicate the existence of an axion like particle, though at the moment it is difficult to reconcile any 
standard axion like particle with these results. For a summary of some of the possible ways out to 
this difficulty refer to [T3] . 

In this note we would like to study some cosmological implications/effects in the context of axions. 
To the best of our knowledge, all the existing studies have ignored the finite temperature effects to 
axion decay constant while discussing issues like relic axion abundance etc. We intend to explore 
the possible impact of such corrections on axion physics at different epochs of the early universe. 
Before proceeding with the finite temperature effects, we fix our notation. It is worthwhile to 
mention that axions share the basic properties with the pions, notably, the axion coupling to two 
gluons is given by the anomaly term (the two photon coupling has an analogous form) 

L-jgcret,, (i) 

where a s is the QCD coupling constant and G b is the dual of the gluon field defined as 

G fiv = ~^ e iivpoG p , (2) 
where b is the color index. Let us rewrite the above as 

L = Qagg^pvpcjG ^ G P O , (3) 

where g agg is the axion coupling which is inversely proportional to /„. The axions can be further 
divided into two categories - (a) Hadronic axions which couple to quarks only and (b) non-Hadronic 
axions having a coupling to the leptons. Mass of hadronic axions is given by 

m ^TTrjr- (4) 

where f n = 93 MeV is the pion decay constant, and z = m u /md is the mass ratio of up and down 
quarks, z can lie within the interval 0.3 — 0.7 [fSJ. Here we take z = 0.56 jTBJE]- 

Having noted that the axion coupling to massless gauge bosons is the same as pion-photon-photon 
coupling, we model the temperature effects in an analogous way as well. To this end, we first note 
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that the results on chiral dynamics and anomaly at finite temperature IT9] imply the following 




(5) 



where f T = f n (T = 0) ~ 93 MeV is the zero temperature pion decay constant measured in the 
laboratory. It is important to keep in mind that such a result is valid only for T < In the case of 
axions, it is natural to look at temperatures lower than the axion decay constant as it is only then 
that the axions can be thermally produced. Therefore, the above result for the pions can be directly 
extended to the axion case. The analysis in [19 however brings out a crucial point. The effective 
pion-photon-photon coupling is inversely proportional to the decay constant and one would have 
expected that the coupling changes the way dictated by the change in decay constant with the 

temperature, i.e., one would have expected the coupling to change as g n77 — > g nll (l — ^j?) or 



that to 0(T 2 1 f%), g^{T) ~ 1/ f n (T). What is found instead is that to this order g nll (T) ~ f-wiT), 
i.e., 



Note that since the limits on axion mass indicate a very small mass, we do not consider temperature 
dependent effects for axion mass. 

We are now ready to explore two epochs in the history of the universe, particularly important 
and relevant for axion cosmology. First, we study axion thermalization and abundance in the very 
early universe including the temperature effects. The second case investigated is the axion hadron 
interactions in the post QCD era. In both these cases we find that the temperature effects - via 
axion or pion decay constant, pion mass or axion-gluon-gluon coupling, whichever is applicable at 
the relevant epoch - do lead to change in the results for certain values of the temperatures. 

2 Axion thermalization and abundance 

As mentioned earlier, axions can be thermally produced in the universe once the temperature falls 
below the axion decay constant. However, axions can also be produced via some non-thermal 
processes (see [23] and references therein). In |24j . the authors identify the conditions under which 
there is significant production and/or thermalization of axions. We therefore closely follow [21] 
but now include temperature effects as mentioned before. The main thermalization processes 
are axion-gluon and axion-quark scatterings. From Eq.(j3J) and Eq.©, we see that both f a and 

1 It is worthwhile to point out that the configuration considered in does not corespond to the physical 
situation while the authors in 20] consider the decay of pseudoscalar to real photons and find difference with the 
results obtained in |19|. See f° r a discussion about this difference and the resolution. However, we only use the 
results concerning the change in decay constant with temperature and therefore this issue does not affect us here. 
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We therefore take this form for the axion-gluon-gluon coupling and work with this. 1 
Further, following ^211221 ; we have for the pion mass (again in the regime T < /„-), 




(7) 
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Qagg decrease with the inclusion of the temperature effects. This means that the axions can be 
thermally produced at a slightly lower temperature and the interaction rate also goes down a bit. 
The Boltzmann equation for the abundance in an expanding universe is given by 

where x = Jfi, H is the Hubble expansion rate and T is the thermally averaged rate of reaction of 
axion processes. In the radiation dominated era, the Friedmann equation yields 

V 45 J M P ' K J 

where g e ff is the effective degrees of freedom at the temperature T. For a generic process a + i — > 
j + k k ^ a), thermally averaged rate is 

r = — / d Pa d Pl d Pj d Pk (2n) 4 s\ Pa + Pi - Pk - Pk )\M\ 2 Ui{\ + fj){\ + f k ) , (io) 

where d P = i 2 t)^2E ' are ^ ne equilibrium distribution functions and n a is the equilibrium axion 
number density. Compared to the expression for the thermal averaged interaction rate in [21], we 
will have an extra temperature dependent factor stemming from the correction to the axion-gluon- 
gluon coupling and the rate now looks 

Y eq is the abundance when axions are in thermal equilibrium with other SM particles. In terms of 
variables 77 = and k = Xjj-, Eq.© becomes 



Retaining only the leading term, the solution is 



dx 



i](x) — 1 + C ex P 



k 



x 18x 3 



(13) 



where C is a constant determined by boundary condition, 77 = for = 1. This gives 77 = 
at x ~ 1.08 instead of at x — 1 as the boundary condition and the constant C is evaluated 
accordingly. FigO] compares the axion abundance with and without the temperature effects taken 
into consideration for f a = 1.2 x 10 12 GeV. We see that there is a perceptible change in the 
abundance for the temperature range (6 — 10)xl0 11 GeV. As a consequence of the shift in the x 
value from unity to 1.08, the results in [21] shift by this factor. 



3 Hadronic axion in the post QCD era 

The hadronic axion has no lowest level coupling to the charged leptons and therefore the induced 
coupling is expected to be small enough to be neglected altogether. For updated bounds on hadronic 
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Figure 1: Abundance of axions for f a = 1.2 x 10 12 GeV. 

axion in this context refer to |2*5] . Although the exact limits on the mass and decay constant of 
the axion are model dependent, up to order unity factors, the limits from various sources suggest: 
f a > 0.6 x 10 9 GeV and m a < 0.01 eV. Stellar energy loss constraints suggest that in the post 
QCD era, the axion coupling to electrons and photons are not significant. The axion coupling to 
nucleons is constrained by the observed neutrino signal from Supernovae 1987A [5|l29j. Supernovae 
observations give stringent constraints on the mass and coupling: m a > 0.01 eV or f a < 0.6 x 10 9 
GeV. However it still leaves a domain for hadronic axion where f a lies in the range 3 x 10 5 GeV 
and 3 x 10 6 GeV |2*6*1 |2"T]. But this range is disfavored by a detailed combined analysis of hadronic 
axion with cosmological data like large scale structure, cosmic microwave background, supernova 
luminosity distances, Lyman-a forest and Hubble parameter , which obtains the following limits: 
m a < 1.05 eV or equivalently f a > 5.7 x 10 6 GeV. However, in the post QCD era of interest to 
us here, the number density of the nucleons falls much faster than that of pions and therefore, 
nucleons can be effectively neglected compared to the pions. That the nucleon contribution to the 
thermalization and decoupling of axions in this era is much smaller than the pion contribution has 
been explicitly verified in [2*o] . 

The axion-pion interaction is of the form [26j 

L a7T = yj- (^Tr+^TT- + n°n-d,n + - 27r+7r-^7r ) d^a . (14) 

Jit J a 

In hadronic axion models, the coupling constant is 

CaK = WVz) • (15) 

The relevant processes are an^ — > tt ^ and an — > tt + 7t~ . Following [213 ESj - we calculate the 
axion decoupling temperature in this context but now taking the finite temperature effects in 
the pion decay constant and mass into account given in Eq.(j5|) and Eq.((7j) respectively (valid 
for T < f n ). The decoupling temperature is that where the Hubble expansion rate equals the 
thermally averaged interaction rate. In the Fig|21 we plot the expansion rate of the universe (H) 
and the interaction rates (r) - with and without temperature effects put in. The points where 
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the T curves cut the H curve give the decoupling temperature. This is shown in the figure for 
f a = 10 7 GeV. It is clear from the figure that the decoupling temperature is lowered once the 
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Figure 2: Hadronic axion reaction rate with(solid) and without (dashed) temperature effect for 
f a = WGeV. 

finite temperature effects are included and the reason is easy to see. In the epoch of interest, the 
temperatures are much lower than the axion decay constant and therefore, we do not consider any 
temperature dependent effects in f a , while f w decreases and increases with temperature. Both 
these effects lead to larger interaction rate and therefore, the axions decouple later in time or at 
lower temperatures. From the figure, it is also clear that the difference in the two cases is not 
very large. We also compare decoupling temperature of axions in the present calculation with the 
previous ones (without the temperature effects) for several representative values of f a in Table ^ 
We find that for higher values of f a , the temperature effects start becoming significant. However, 
within the approximation adopted here, it is not possible to explore the range f a > 1.5 x 10 7 GeV 
as the temperatures exceed the pion decay constant, thus taking us away from the validity of the 
approximation employed. It is therefore important to investigate the temperature effects going 
beyond this approximation and exploring the consequences. 



/a(GeV) 


Tdi (MeV) 


T D2 (MeV) 


3 xlO 5 
1 xlO 6 
3 xlO 6 
1 xlO 7 

1.2 xlO 7 

1.3 xlO 7 


26.43 
35.34 
49.84 
81.04 
87.61 
90.1 


26.43 

35.34 

49.5 

79.12 

85.48 

87.9 



Table 1: Decoupling temperature of axions Tdi (without) and Td2 (with) temperature effects for 
different values of f a . 
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4 Discussion 



We have investigated the impact of temperature effects, though in a limited sense, on the axion 
cosmology. We have focused our attention of the temperature effects in the axion decay constant, 
pion decay constant and pion mass in two important epochs of the universe. The results show that 
there can be perceptible differences in the predictions of abundances and decoupling temperatures. 
However, we must remark that these effects have been considered in the approximation when the 
temperature of the universe is lower than the decay constant(s). In the case of axion thermalization 
in the early universe this has a physical meaning because it is only then that the axions can 
be thermally produced. In this case we see some change in the abundance of the axions with 
temperature. In the case of axion interacting with pions in the post QCD era, we find that 
temperature effects in the pion mass and decay constant lead to a lower decoupling temperature 
for the axions. We find that as the axion decay constant is increased the difference between the zero 
temperature and finite temperature cases starts becoming significant. However, we have worked in 
the limit when temperatures are smaller than /„.. This puts a natural limit to the region that we 
can explore. It is worthwhile to try to explore what happens when T > f n . This will require going 
beyond the validity of the present calculation. We would also like to mention that a more detailed 
and consistent calculation requires computing the matrix elements including thermal loops also to 
the same order in temperature corrections and then studying the full impact. 

In summary, we can say that we have studied the impact of temperature effects due to a very limited 
source (s) entering the axion cosmology. The results show perceptible departure from the zero 
temperature calculations. Let us just compare the situation with the case of neutrino decoupling. 
In that case, naively, the temperature effects shift the decoupling temperature by a seemingly small 
amount. However, for precision studies, such a shift makes a significant difference and hence turns 
out to be important. Our approach is also guided by this fact. See [2U| (and references therein) for 
related discussion. Whether, finite temperature effects are actually significant and substantial can 
be decided by more detailed computations, taking into account all sources of corrections to this 
order. Also, in the case of axion decoupling in the post QCD era, it'll be interesting to try to go 
beyond the approximation of T < f m and include finite temperature effects in the masses and all 
other parameters entering the expressions. This can have serious implications for axion cosmology 
and can lead to very different bounds/limits on various parameters, thereby possibly changing our 
current understanding. 
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